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riiis  report  presents  the  derivation  of  an  equation  lor  the  radar  haikscatter  coellicieni  Iroin 
a tilted,  slightly  rough  surface.  The  basic  technique  employed  in  the  solution  is  the  small 
perturbation  inetliod.  along  with  the  Fourier  Iransform.  The  incident  wave  is  allowed  to  have 
an  arbitrary  polari/ation.  and  specific  results  are  shown  for  hori/ontal.  vertical,  circular,  and 
elliptical  polarizations.  One  conclusion  of  this  report  is  that  the  innuence  of  a given  slojx- ol 
the  tilted  rough  surface,  upon  the  calculation  of  the  radar  backscatter  coefficient  depends 
upon  the  pol.iri/ation  of  the  incident  electromagnetic  wave. 
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SUMMARY 


The  work  described  in  (his  report  was  done  under  two  separate  but  related  tasks. 
One  task  is  that  of  generating  radar  image  simulations  of  various  types  of  terrain.  To 
adequately  simulate  a radar  image,  appropriate  radar  scattering  models  are  needed  to 
determine  the  final  image  gray  tone.  The  particular  model  developed  in  this  report  is 
that  of  a tilted,  slightly  rough  surface.  A surface  such  as  this  could  simulate  a road, 
runway,  desert,  or  any  other  surface  that  meets  the  constraints  of  the  solution.  The 
second  task  is  obtaining  military  geographic  information  using  radar.  If  a calibrated 
radar  were  used  to  provide  appropriate  radar  measurements  of  terrain  and  if  these 
measurements  were  coupled  with  an  applicable  scattering  theory,  then  physical  and 
electromagnetic  parameters  that  characterize  the  surface  could  be  calculated. 

Therefore,  this  report  presents  the  necessary  mathematics  to  solve  the  problem  of 
electromagnetic  wave  scattering  from  a tilted,  slightly  rough  surface.  The  polarization 
of  the  incident  wave  is  allowed  to  be  arbitrary,  and  specific  results  are  shown  for 
horizontal,  vertical,  circular,  and  elliptical  polarizations. 
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BACKSt  ATTtRING  OF  RADAR  WAVtS  FROM 


A TILTED.  SLIGHTLY  ROUGH  SURFACE 

INTRODUCTION 

Purpose.  This  report  develops  an  equation  for  the  radar  hackscatter  eociTicicnt 
from  a tilted,  slightly  rough  surface.  I he  basic  technique  employed  in  the  solution  is 
the  small  perturbation  method,  along  with  the  Fourier  transfomi. 

Background.  The  problem  of  treating  radar  wave  scattering  from  random  rough 
surfaces  has  been  studied  for  many  years.  There  are  many  practical  applications  for  a 
theory  which  could  explain  this  problem.  One  application  is  the  development  of  the 
basic  understanding  of  scattering  phenomena.  Another  application  uses  scattcnng 
theories,  along  with  certain  radar  measurements,  to  detemiine  quantitatively  various 
physical  parameters  of  the  surface  itself.  For  instance,  if  one  has  ex|)erimentally  mea- 
sured the  radar  buckscatter  coefficient  at  several  polarizations  or  frequencies  over  a 
particular  terrain  surface,  one  may,  by  using  an  appropriate  theory,  he  able  to  calculate 
surface  roughness  properties,  moisture  content,  or  other  surface  characteristics.  Still 
another  application  is  the  area  of  radar  image  simulations.  If  one  is  interested  in 
making  a radar  image  simulation  of  a particular  feature,  then  theoretical  radar 
scattering  models  can  be  used,  along  with  the  radar  range  equation,  to  obtain  quantita- 
tive expressions  for  the  simulated  radar  image  gray  tones.  The  above-mentioned  appli- 
cations provide  incentive  for  finding  solutions  to  various  radar  scattering  problems. 

A problem  in  radar  scattering  that  needs  theoretical  development  is  that  of 
scattering  from  a tilted,  slightly  rough  surface.  One  solution  to  this  problem  has  been 
developed  by  Valenzuela.'  However,  his  final  result  considers  only  the  effect  of  a 
surface  that  is  tilted  in  the  direction  orthogonal  to  the  plan*  of  incidence.  In  general, 
there  will  be  two  components  of  slope  that  must  be  considered,  one  component  of 
slope  in  the  plane  of  incidence  and  the  other  component  in  a direction  orthogonal  to 
the  plane  of  incidence.  The  following  derivation  will  consider  the  scattering  of  electro- 
magnetic waves  from  a slightly  rough  surface  that  has  an  arbitrary  tilt.  The  derivation 
will  be  confined  to  a consideration  of  waves  in  the  backscatter  direction  only;  however, 
the  incident  polarization  will  be  taken  as  arbitrary.  The  final  result  of  the  derivation 
will  be  an  equation  for  the  radar  backscatter  coefficient  in  terms  of  the  two  slopes. 
The  small  perturbation  technique  will  be  used,  and  results  up  to  first  order  will  be 
considered.  The  derivation  uses  the  rationalized  MKS  system  of  units.  Also,  results 
published  previously  will  be  used  to  obtain  a final  solution.’ 

G.  R.  Valenzuela.  "Scallerinp  of  I'lectromapnetic  Waves  from  a Tilled.  Slightly  Rough  Surface."  Radin  Sciencr. 
^ Vol.  3,  No.  1 1 , November  1968. 

R.  A.  tlevenor,  Backscaucriw  of  l.lectromagnctic  Wares  from  a Surface  Composed  of  Two  Types  of  Surface 

Rounhness,  KTL-rR-71-4,  U.S.  Army  I ngineci  Topographic  Lahoialorics.  l-ort  Ucivoir,  VA..  Uclobcr  1971. 

\ 11-737  675. 
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ANALYSIS 


The  geometry  of  the  problem  to  be  studied  in  this  report  is  given  in  f igure  I . 
A right  handed  rectangular  cartesian  coordinate  system  (x.  y,  /)  is  set  up  with  the  x/ 
plane  coinciding  with  the  plane  of  the  paper.  The  y axis  is  perpendicular  to  the  plane 
of  the  paper  and  pointing  into  it. 


A plane  wave  with  a time  harmonic  of  expfju^t)  is  incident  from  free  space  onto  a 
dielectric  surface  S.  The  xz  plane  will  be  taken  as  the  plane  of  incidence  with  Di  a 
unit  vector  in  the  direction  of  propagation.  An  equation  for  Q|  can  be  written  as 
follows; 


Oi  = sin  0 - cos  6 ( I ) 

where  4^  is  a unit  vector  in  the  x direction. 

a^  is  a unit  vector  in  the  z direction. 

0 is  the  angle  of  incidence  measured  with  respect  to  the  z axis. 

The  line  underneath  a symbol  shall  be  used  to  represent  a vector  quantity.  The 
surface  S will  be  considered  equal  to  the  sum  of  two  surfaces;  one  of  which  is  a plane 
and  the  other  of  which  is  a random,  slightly  rough  surface.  If  wc  represent  the  z co- 
ordinate of  the  surface  S by  p (x,  y),  then  we  can  write  the  following  expression  for 
P (X,  y); 
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p (X.  y)  = Z (X. y) + s(x.  y) 


(2> 

where  Z (x,  y)  represents  an  expression  for  the  l coordinate  of  a plane  that  has  an 
arbitrary  tilt,  and  s(x,  y)  represents  an  expression  for  the  z coordinate  of  a random, 
slightly  rough  surface. 

The  surface  function  s(x,  y)  will  be  considered  as  being  generated  by  a wide 
sense  stationary  gaussian  random  process  with  a zero  mean,  a variance  o^,  and  a 
correlation  function  0 (t  , r ).  The  unit  vector  Dj  is  in  the  direction  of  the  receiver. 
For  the  case  of  a monostatic  radar  system.  D:  is  simply  equal  to  the  negative  of  0|. 
The  unit  vector  n is  an  outward  normal  to  the  plane,  Z (x,  y)  and  can  be  written  as 


.h. 

VTTzfTzJ” 

" 7 


(3) 


where  Z,  = — and  represents  the  slope  of  the  plane  in  the  x direction. 

' ox 

3Z 

Z = and  represents  the  slope  of  the  plane  in  the  y direction, 
y oy 

a^ , Qy , and  are  unit  vectors  in  the  x,  y,  and  z directions,  respectively. 

The  local  angle  of  incidence,  O',  is  defined  as  the  angle  between  n and  -q,  and  can 
be  calculated  as  follows: 


cos  0'  = - n,  • n 


The  incident  electric  field  (Ej)  can  be  written  as 

c _ -jkoOi-i 

Ei  = e 


(4) 


(5) 


where  e^  is  an  arbitrary  unit  polarization  vector  and,  in  general,  is  complex. 


is  the  free  space  propagation  constant, 
r is  a position  vector  equal  to  xa  + ya  + za  . 

X y 


The  polarization  vector  Cj  is  chosen  as  arbitrary  to  consider  the  results  for  a number  of 
different  incident  polarizations,  such  as  horizontal,  vertical,  and  circular.  The  time 
harmonic  portion  of  the  wave  has  been  suppressed  and  will  not  be  needed  in  the  rest 
of  the  analysis. 
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A locul  rcctunt;ular  coordinate  system  (x'y'z't  is  set  up  such  that  the  z‘ axis  coin- 
cides with  the  unit  normal  d-  I he  y'  axis  (into  the  pa|>er)  is  placed  such  that  it  is  pet 
pendicular  to  the  local  plane  of  incidence  formed  by  n and  ni  ■ An  equation  for  a unit 
vector  in  the  y'  direction  (a'  ) would  he  as  follows; 


= 0 X D|/(lo  X u,|) 


(6) 


When  the  definitions  for  o,  and  o as  given  by  equations  ( I ) and  (3)  are  placed  in  (6t, 
the  following  expression  for  a’  results: 


i>'  = la^  Zy  cosO  + a^  (sinO  - cos0)  + a,  Z^  sin0|  D,  (7) 


where 


D,  = I Z^y  + (sin  0 - cos0)^  | 

Since  the  direction  of  the  /'  axis  coincides  with  q.  a unit  vector  (aj ) in  the  direction  ol 
z'  will  simply  be  equal  to  q,  and  a unit  vector  in  x'  can  now  be  calculated  easily. 


ilx  = “y  * ifi 


al  = / a,  I sin  0 - Z,  COS0  + Z^sin01  - a„  Z.,  |cos0  + Z,  sin0| 
xiA  X y y>  X 

+ ilj  (Z^  (sin  0 - COS0)  - Z*  cos0|  | Do  D, 


(8) 


(yi 


where 

'>■>-{'*  * 'I  \ 

The  problem  of  electromagnetic  wave  scattering  from  a slightly  rough  surface 
with  no  tilt  has  been  developed  for  both  horizontally  and  vertically  polarized  incident 
waves.  If  the  electric  field  incident  onto  the  tilted  slightly  rough  surface  is  broken  up 
into  local  horizontal  and  vertically  polarized  components,  then  existing  results  can  be 
used  to  obtain  a solution.  The  local  horizontally  polarized  component  can  be  found 
by  computing  the  dot  product  of  with  3^. 


fc-i  • ay  = ay  • c 


■Ao".  •' 


< 10) 


In  general,  the  polarization  vector  g.  is  made  up  of  three  orthogonal  components. 

e.  = e.  a + e.  a e a (II) 

*4  IX  “x  ly  -y  it.  -7 
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I he  local  hon/ontally  polan/cd  cuiiipunent  of  the  incident  field  then  becomes 

t,  (i.’i 

where 

a,  = |Zy  (e,^  cost)  + e,^  sinO)  + e,^  (sinO  - /.^  cosO)|  I), 

Kquatiun  (12)  now  provides  an  expression  for  the  local  hori/.ontally  polari/.ed  incident 
electric  field  in  tenns  of  the  two  slopes  of  the  plane,  the  angle  of  incidence  and  the 
components  of  the  polarization  vector. 

To  compute  the  local  vertically  polarized  component,  we  will  work  with  the 
magnetic  field.  The  incident  magnetic  field  (U,)  can  be  calculated  as  follows: 

n»  X t, 

M,  = (13) 

where  a;  = 2rrf  and  f is  frequency. 

= magnetic  permeability  of  free  space. 

The  local  vertically  polarized  component  of  the  incident  wave  can  be  determined  by 
computing  Hj  • a^. 

Hj  • a'  = aj  e ' (14) 

where 


^2  “ [Z  Cjy  - (sin  0 - Zjj  cos  0)  (ejj  sin  0 + e^j^  cos  0)1  D, 


The  components  of  the  total  electric  field  in  the  free  space  medium  for  the  case  of  the 
local  horizontally  polarized  incident  wave  were  taken  from  ETL-7R-71-4  and  are  pro- 
vided below; 

E^-  = y"  y*  A^i  (k^,  ky)  exp  (jk^x' + jk^y'  - jk^  z' ) dk^  dky  (15) 


^ R.  A.  Kevenor,  Backscattermg  of  Electromagnetic  Waves  from  a Surface  Composed  of  Two  Types  of  Surface 
Roughness.  ETL-TR-71-4,  U.S.  Army  Engineer  Topographic  Laboratories,  I'ort  Belvoir,  VA,  October  1971, 
AP-7J7-67S. 


5 


a,  exp  (-jk^x  sin  0 ) je 


+ R,e 


•f 


(k^.ky)exp(jk^  x’+jky  y'-jk^  /)dk^  dk. 


(k^,ky)exp(jk^x'  + jkyy'  - jk^z')  dk^  dk^ 


k,  = \/<  - kl  - k’, 

k cos  o'  - \/k'^  - k^  sin*  O' 

i! 

kjj  cos  o'  + k'^  - k^  sin^0' 

K\  <kx-ky)  = a,  k^ky  O' 

Ayi  (k^.ky)  = - a,  (k^  + k^  k^)  O’ 

\\  = - a.  ky  k;  O' 


jT,  (k'^  -kl)  S (k^  + k„  sin  O',  ky) 


S(k^.k^) 


- > f f i ' ■’‘‘x*'  J - I < 

= — f j s(x,y)c  * c ’ dy  dx 


The  det'inition  of  the  various  parameters  given  above  are  provided  in  the  back  of  this 
report.  Tlie  terms  with  the  double  integrals  in  equations  ( 1 5)  through  (17)  represent 
the  scattered  fields,  however,  the  tenn  in  equation  ( 16)  without  double  integrals  repre- 
sents an  incident  wave  plus  a specular  component.  The  parameter  s(x',  y')  is  the 
expression  for  the  slightly  rough  surface  in  the  primed  or  local  coordinate  system. 


The  components  of  the  total  electric  field  in  the  free  space  medium  for  the  case 
of  the  local  vertically  polarized  incident  wave  are  taken  from  hTL-TR-71-4  and  given 
as  follows;*' 


0 

E,.  . a,  n.-*" cos «•  { R„  e-*”'  "'*' 


If 


\2  ‘**‘y 


(18) 


Ey'  = 


If  \2  ‘^'‘y 


(19) 


n _ -jk  x'sinO'  . , ( jk  z' 

E^.  = - Tj  e ° sm  0 ) e ° 


cose'  „ -jk„z’cose 

+ R„  e ° 


} 


If 


\2  ‘^'^y 


(20) 


where 


R„  = 


k'^  cos  6'  - k„  V k'^  - k^  sin^  6' 


k'^  cos0'  + k„Vk'^  - k^  sin^  6' 


R.  A.  Hevenor,  Backscattering  of  Electromagnetic  Waves  from  a Surface  Composed  of  Two  Types  of  Surface 
Roughness,  ETL-TR-71-4,  U.S.  Army  Engineer  Topographic  Laboratories,  Fort  Belvoii,  VA,  October  1971, 
AD-737-675. 


J 
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= ' { ^1^0  Qi  k^^  Q,  +k^  k^'  k^  k'^W  + k^  k'  k^  k'^  vj 

®o  = K ^'z  K - >7  l<o  K K Q> 

a,.  = V kj  k;  k'^  + n kl  k;  k'^  ^vKKK^rikl  k'/k^ 

a, 2 = 7J  ky  k'^  k'"  k^  + T)kl  ky  k^  k^ 
a^,  = k'  ky  + 7?k;,  k^  k^  ky 

^ R.  A.  Hevenor,  Backscattering  of  Electromagnetic  Waves  from  a Surface  Composed  of  Two  Types  of  Surface 
Roughness,  UTL-TR-71-4,  IJ.S.  Army  I nginecr  Topographic  Laboratories,  Lort  Belvoir,  VA,  October  1971 
AD-737-675. 
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The  parameter  R in  equation  (21)  is  simply  the  distance  from  the  origin  of  the  local 
coordinate  system  (x'y’z')  to  the  receiver  in  the  far  field,  liquations  (22)  through  (24) 
can  easily  be  evaluated  using  previous  results. 


(ky  sin  d',  0)  = 0 


A^2  0)  = aj  T?  A cos  6'  S (2  sin  O',  0) 


'.x2 

where 


(25) 

(26) 


A ' - { k.  k'  k?  c„  k k,  k;  k„  k'»  c,  + k,  k'.  k„  k’=  / (k;  k;  k'‘  k k’  k';  k. . 

= j(ko  -k'")(l  -R||)cos07(27rkJ 
C,  = j (2k^,  sin^  0')  T||  (k^  - k'^  )/(27r  k'^ ) 

Cz  =-jT||  ky  sin^  0' (k^  - k'^ )/(27r  k'^ ) 

Ay,  (k^  sin0',O)  = j a,  Tj^(k^  -k'^)S(2k^  sin  0',  0)/ 1 27r  (k^  + k'j,)I 


Ay2  (k^  sin  O',  0)  = 0 


A^j  (k^  sin  O',  0)  = 0 


A^^  (ko  sin  O',  0)  = aj  tj  A sin  O'  S (2k^  sin  O',  0) 

The  backscattered  far  field  can  now  be  written  as  follows: 

-ik^R 


(27) 

(28) 

(29) 

(30) 


Ej  = (27r)^  2 cos  O' 


(jk„)  e 


47rR 


I ^22  ilz  ^33^ 


S(2k^  sin0’,  0)(31) 


where 


A,i  = a2  1?  A cos  O' 

A32  = ja,  T^(k^-k'')/[27r(k^+k;)] 

A33  = 32  77  A sin  O' 

By  using  previously  derived  results,  one  can  now  write  in  terms  of  the  unit  vectors  of 
the  original  reference  coordinate  system  (x,  y,  /.). 
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-jk^R 

= Ctt)^  2 cosO'  I + Ay  Uy  + S(2k„  sinO',0)  (32) 


K ^ “xl  * “yl  ■*■  “z 


Ay  = A,,  a^2  + A22  Oy,  + A33 


\ " '^"  “x3  “y3  -^33  “z3 


^x  = “xl  ^x  +“x2  ^y  + “x3 


Oj^i  = (sin  6 -Z^  cos  (?  + Zy  sin  6)  D, 

C(^2  ^ ■ ^y  ® ■‘'  ^x 

a , = [Z  (sin  0 - Zj^  cos  0)  - Zy  cos  0|  Djj  D, 


Uy  = ay,  a^  + «y2  ay  + «y3 


tty,  = Zy  D,  COS0 

ttyj  = (sin  0 - Zj^  COS  0)  Di 
“y3  ^ ^y  ^ 


= a^,  + a^^  Sy  + “^3 


“zi  = - Z,  Do 


tt  -X  — - Z D 

z2  y o 


«z3  = 


If  the  radar  receiver  is  sensitive  to  an  electric  field,  which  has  a unit  polarization  vector 
e^,  then  the  received  field  (E,^ ) will  be 


F = e • E 


The  unit  polarization  vector  e^  will  be  allowed  to  have  complex  components  in  all 
three  dimensions. 


C = C 11  + C 11  + C 11 

rx  -X  ry  -y  n ~i 


The  received  electric  field  can  now  be  written  as 


-jk  R 

(jkJe  “ 

= 2 (2?:)^  cos  e'  - - App  S (2k^  sin  O',  0) 


where 


App  = e^^  + Ay  e^y  + A,^  e,^ 


The  radar  backscatter  coefficient  0“  can  be  computed  from  the  received  electric  field 
as  follows: 


O = ^r> 


The  quantity  A^^  is  the  area  of  the  surface  S which  is  illuminated  by  the  radar  beam. 
The  brackets  around  E^E*  are  used  to  indicate  the  calculation  of  a statistical  average. 
The  random  portion  of  E^  comes  from  S(2k^sin0',  0)  since  the  surface  function 
s (x',  y')  is  random. 

k^  cos^  s' 

<E  E*>  = A„  A*  <S(2k  sinO',0)S*(2k  sinfl',0)> 

It  remains  now  to  calculate  the  statistical  average  indicated  in  the  above  equation. 


<S(k^,ky)S*(k,,ky)>  = <S(k,,ky)S(-k^,-ky)>  = 


If  If 


-jk  (x'-x")  -jk  (y'-y") 

e dy  dx  dy  dx 


<s(x',y')s(x",y")>  e 


If  we  let  T = X - X and  t = y - y and  we  realize  that  <s  (x  , y ) s (x  , y ) > is  the 

X y 

correlation  function  of  the  surface,  then  we  have 
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<s(k^.k,isi-k^.-k,i>  =.  ydx'y  dy-y  d,,y  d-, 


0 (r^ . r ) c ' 0 ' 


The  tntei’rals  in  x'  and  y'  at  first  appear  to  be  meaninj’less.  however,  the  result  of  the 

two  integrals  ean  be  represented  approximately  by  the  illuminated  area  , as  long  as 

the  dimensions  of  this  area  are  much  greater  than  the  correlation  distances  associated 

with  0 (t^,  Ty ).  The  other  two  integrals  m and  can  be  used  to  define  the  surface 

roughness  spectrum,  W (k  , k ): 

^ y 


W(k^,ky 


where 


/ 


dTy  C ) e e 


0(T^.ry)  = o,  C(r^,T^) 

The  result  fpr  the  needed  average  can  now  be  written  as 


<S(k  •,  k )S*(k  ,k  )>  = 


A„ 


W(k  ,k  ) 

X ’ y 


The  radar  backscatter  coefficient  can  now  be  written  in  its  final  form: 
a'^  = Hrr^  k^  cos^  0 A„_  A*  W (2k  sin  o\  0) 

O I PP  PP  O 


This  answer  is  a function  of  two  slopes  and  ) and  the  surface  roughness  spec- 
trum. In  the  next  section,  we  will  pick  a specific  form  for  the  correlation  function, 
and  the  roughness  spectrum  will  be  computed.  This  in  turn  will  allow  us  to  calculate 
a°  as  a function  of  the  incidence  angle,  0.  Before  we  do  this,  however,  two  special 
cases  of  slope  must  be  considered  separately  so  that  computational  difficulties  will  not 
arise.  The  first  special  case  occurs  when  Z = 0 and  Z = 0.  It  can  be  seen  from  the 

I o ^ 

equation  for  a^  (equation  7)  that  at  0 = 0 we  will  encounter  a zero  over  zero  condi- 
tion when  both  slopes  vanish.  Therefore,  for  this  case  we  will  need  to  calculate  new 
values  for  a, , a„ , a,,  a, , and  aj  ■ In  this  special  case,  the  unit  vectors  a,  and  n are  both 

equal  to  a,.  The  quantity  a then  becomes 
^ y 

I 

a = a 
-y  ‘V 
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This  makes  a equal  to  a . The  parameters  ai  and  a^  tan  be  calculated  easily: 


i.  = a;  • { 


u + c ii  + c a 
ly  “y  12  “2 


aj  = - (e  sin  0 + e.  cost?) 

The  second  special  case  occurs  when  = 0 and  7.^  = tan  0.  It  can  be  seen  from  equa- 
tion 7 that  this  condition  again  yields  an  indeterminant  zero  over  zero  situation.  We 
will  define  a unit  vector  Pa  that  will  be  used  to  calculate  a ^ ; 

n3  = a^  ‘-'OS  0 + sin  6 

The  quantity  a’y  will  now  be  defined  in  the  following  expression  for  this  special  case: 


a'  = g X 03 

' ^ I n X n,  I 


n X D3  = 


a = a 
-y  -y 

Since  is  equivalent  to  n,  we  can  now  calculate  for  this  special  case. 


3x  = 


3x  3z  ® 


V 1 + tan^  e 


The  only  remaining  computation  that  needs  to  be  performed  is  to  detennine  a,  and  aj ; 


=ay  • {‘-'ix  a, 


+ C-  a + C-  a 
X >y  ^ z 


a,  = - 


(Cj^  sin  6 + cos  0) 


It  can  now  be  seen  that  the  solutions  for  a,  and  for  this  special  case  arc  identical 
to  the  solutions  for  aj  and  aj  for  the  first  special  case  given  previously,  which  was  not 
obvious.  This  completes  the  derivation  of  the  radar  backscatter  coefficient  as  a func- 
tion of  the  two  slopes  Z and  Z . The  final  result  does  not  consider  depolarization 
due  to  the  slightly  rough  surface  itself,  since  only  first  order  terms  in  perturbation  were 
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used.  However,  depolarization  due  to  the  tilted  plane  is  considered.  In  the  next 
section,  sample  calculations  will  be  made  for  the  specific  case  where  C t^)  has  a 
gaussian  form. 


RHSULTS 

The  purpose  of  this  section  is  to  evaluate  the  radar  backscatter  coefficient  (equa- 
tion 36)  for  a particular  surface  correlation  function.  There  are  two  correlation  func- 
tions that  are  used  most  often  in  theoretical  scattering  work,  and  these  are  the 
exponential  and  the  gaussian  functions.  If  we  use  a gaussian  correlation  function,  then 
C Ty ) will  have  the  following  form: 


C(T^,Ty) 


-(T 


where  C is  the  correlation  distance. 


(37) 


In  using  the  above  expression  for  C (t  , t ),  we  have  assumed  an  isotropic  correlation 
function.  This  assumption  would  not  be  valid  for  surfaces  that  contain  some  type  of 
periodic  structure  in  one  direction,  i.e.  a plowed  field.  Recalling  the  definition  of  the 
surface  roughness  spectrum  W (k^^,  k^ ) and  inserting  equation  37  into  that  expression 
will  yield 


W(k^,ky)  = 


-a’+k’ )tV4 

e * y 


A computer  program  was  written  for  the  calculation  of  equation  36  using  the  surface 
roughness  spectrum  given  above.  The  computer  program  listing  for  a vertical  polariza- 
tion solution  is  provided  in  the  appendix.  A set  of  numerical  results  for  four  different 
polarizations  is  given  in  figures  2 through  1 3.  Kach  figure  presents  two  graphs  of  radar 
backscatter  coefficient  in  decibels  versus  angle  of  incidence  in  degrees.  The  backscatter 
coefficient  in  decibels  is  related  to  equation  36  as  follows: 

o°  (in  decibels)  = 10  log,o  a° 

The  backscatter  coefficient  on  the  right  side  of  the  above  equation  is  computed  from 
equation  36.  The  numerical  values  used  for  frequency,  correlation  length,  standard 
deviation,  and  relative  dielectric  constant  were  held  fixed  throughout  all  calculations 
and  are  given  on  each  figure. 

The  four  polarizations  considered  were  horizontal,  vertical,  circular,  and  elliptical. 
The  unit  polarization  vectors  for  each  of  the  four  cases  are  given  below: 

for  horizontal  polarization  gj  = = Uy 


\5 


For  vertical  pulari/.ution 
For  circular  polarization 

For  circular  polarization 

For  elliptical  polarization 

For  elliptical  polarization 


Si  = e,  = - (a^  cos  0 + a,  sin  0) 


I ay  - C cos  0 + a^  sin  0)| 

1 ( J"/2  „ I 

Cj  = jay  - C (a^  cos  0 + Uj  sin  0)^ 

^fy-e  (a,  cos 0 + a,  sin  0)^ 

1 ( j»/4  ( 

I ® sin  0)  I 


There  are  three  figures  associated  with  each  polarization.  One  figure  shows  the  effect 
of  changing  only  the  slope  in  the  plane  of  incidence  (2.^).  The  next  figure  shows  the 
effect  of  changing  only  the  slope  in  the  plane  orthogonal  to  the  plane  of  incidence 
(Zy).  The  third  figure  shows  the  effect  of  having  an  equal  slope  in  both  the  x and  y 
directions.  For  all  figures,  a curve  is  plotted  for  the  case  where  no  slopes  exist  (Z^  = 0, 
Zy  = 0).  The  results  for  each  polarization  will  be  discussed  separately  below. 

The  horizontal  polarization  results  are  presented  in  figures  2 through  4.  By  com- 
paring figures  2 and  3,  it  can  be  seen  that  infiuences  horizontal  polarization  much 
more  than  Zy.  The  value  of  slope,  which  is  used  in  all  calculations,  is  equal  to  1.2.  In 
figure  2,  adding  a large  slope  in  x lowers  the  value  of  o°  over  a range  of  incidence- 
angles  extending  from  0°  to  approximately  25°.  For  angles  of  incidence  greater  than 
25°,  the  o°  curve  for  zero  slope  drops  off  much  faster  than  the  curve  associated  with  a 
slope  in  x.  The  crossover  point  in  figure  2 occurs  at  the  point  where  O'  = 0.  The  cross- 
over point  in  figures  3 and  4 are  more  difficult  to  interpret  since  both  local  vertical 
and  local  horizontal  polarization  components  exist.  Figure  4 shows  the  result  of  com- 
puting the  backscatter  coefficient  for  a surface  that  has  a large  slope  in  both  x and  y. 

The  vertical  polarization  results  are  given  in  figures  5 through  7.  By  comparing 
figures  5 and  6,  one  sees  that  Z has  a much  greater  inlluence  on  o°  than  does  Z^. 
This  influence,  however,  is  almost  a bias  as  we  see  that  the  two  curves  in  figure  b have 
essentially  the  same  shape.  Figure  7 presents  the  results  for  a surface  that  has  slopes  in 
both  x and  y.  In  this  figure,  we  see  that  making  = 1.2  and  Zy  = 1.2  lowers  the 
curve  over  the  range  of  incidence  anglesextending  from  0°  to  approximately  63°. 


Tlie  results  of  tlie  circular  polarization  calculations  are  presented  in  figures  8 
through  10.  F'igure  8 shows  the  inlluence  of  a large  slope  in  x and  no  slope  in  y.  It 
can  be  seen  that  again  we  have  a crossover  p^'int  at  approximately  0 = 25°.  Figure  6 
shows  tlie  inlluence  of  a large  slope  in  y and  no  slope  in  x.  I he  a°  curve  for  = 1.2 
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Figure  2.  Study  of  Variations 
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Figure  3.  Study  of  Variation; 
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BACKSCAHER  COEFFICIENT  IN  db 


Figure  5.  Study  of  Variations 
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Figure  6.  Study  of  Variations 


21 


BACKSCAHER  COEFFICIENT  IN  db 


BACKSCAHER  COEFFICIENT  IN  db 


1 


r 

i 

f 


ANGLE  OF  IN 

Figure  9.  Study 


f 


2< 


h 


ANGLE  OF  INCIDENCE  IN  DEGREES 


Figure  10.  Study  of  Z,  and  Z„  Variations 
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Figure  1 1 . Study  of  Variations 
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Figure  12.  Study  of  Variations 
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Figure  13.  Study  of  Z and  Z Variations 
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and  = 0 is  seen  to  be  niuclt  lower  than  the  curve  associated  with  no  slope,  althouyh 
it  does  have  approximately  the  same  shape.  The  eliect  ot  Z^  given  in  l igure  is 
similar  to  the  elTect  exhibited  in  figure  b tor  vertical  polarization.  In  figure  10.  the 
result  of  having  a large  slope  in  both  the  x and  the  y directions  is  seen.  Now,  we  see 
that  the  crossover  point  is  located  at  0 = 40°. 

In  figures  II  through  13,  the  results  for  one  type  of  elliptical  polarization  arc- 
seen.  Figure  1 1 shows  the  intluenee  of  a large  slope  in  x and  no  slope  in  y.  Once  again 
we  have  a crossover  point  at  0 = 2.‘5°.  Figure  1 2 shows  the  effect  of  a slope  in  y and  no 
slope  in  x.  We  can  see  that  even  though  there  is  no  crossover  point  the  two  curves  do 
not  have  the  same  shape.  Figure  13  presents  the  result  of  having  a large  slope  in  both 
X and  y.  In  this  figure,  the  crossover  point  now  occurs  at  approximately  0 = 42°. 

In  this  section,  we  have  considered  the  intluenee  of  surface  slopes  on  the  calcula- 
tions of  the  radar  backscatter  coefficient  for  four  different  polarizations.  It  was  clearly 
evident  that  influences  the  calculation  of  o°  in  a different  manner  than  Z^.  Also, 
we  were  able  to  see  that  the  various  effects  of  slope  cl  anged  from  one  polarization  to 
another.  Since  the  influence  of  slope  changes  with  polarization  and  also  with  the  direc- 
tion of  the  slope,  any  simple  hueristic  approach  probably  could  not  be  developed  to 
adequately  explain  the  effect  of  slope  on  radar  scattering  in  a correct  manner. 

CONCLUSIONS 

The  following  conclusions  can  be  made  from  the  analysis  and  the  curves  presented 
in  the  last  section. 

1 . A theory  has  been  developed  to  explain  radar  backscattering  from  a slightly 
rough  surface  that  has  an  arbitrary  tilt.  The  effect  of  depolarization  due  to  the  tilted 
plane  was  considered,  but  depolarization  due  to  the  slightly  rough  surface  itself  was 
not  considered. 

2.  The  slope  in  the  plane  of  incidence  (Z^^)  influences  the  calculation  of  the 
radar  backscatter  coefficient  in  an  entirely  different  manner  than  a slope  in  a plane 
orthogonal  to  the  plane  of  incidence  (Z^)  for  a given  polarization. 

3.  The  influence  of  a given  slope  of  the  tilted  rough  surface  upon  the  calcula- 
tion of  the  radar  backscatter  coefficient  depends  upon  the  polarization  of  the  incident 
electromagnetic  wave. 

4.  The  results  of  this  mathematical  model  derived  for  electromagnetic  wave- 
scattering  from  a tilted,  slightly  rough  surface  indicates  that  significant  error  could 
occur  if  only  a single  slope  mathematical  solution  is  applied  to  radar  image  simulation 
problems. 


29 


r 


APPENDIX 

A COMPUTER  PROGRAM  LISTING  FOR  THE  CALCULATION 
OF  THE  RADAR  BACKSCATTER  COEFFICIENT 


i 

Sic*isn  ( I Uf'UI  .OUTPUT,  TAPEti=INPult  TAPE6  = 0UTPUT» 

C this  PkUGPAI  COIPUTES  the  «A0AR  nACKSCATfCK  CO£FFir,:ENT 

C F0<  A GLIGhT.Y  <JUGM  TIUTtu  SUkFAoE 

COMPLEX  .4j,QEIX,QEI^,UAl,UA21,QA2,aC0,QCl, 

IQCZ.QAll.OAZ?.  0A3  3.QA44,QA,aAAll.QAA22,QA/.i3  , 

2QAXH,UAYH,QA2H,QALP, QALPC. J3IG0A,QE«X,QC?7 
RFA0(';  2AP,2Y?,XKOL,XKOSIG,eR 

101  FORMAT <6P13.1 ) 

PI=3. 1415926535898 

U»4.0E-7 

U0=PI»U 

E0=8.e54F-12 

XK=2.*PI»F ’SORTIUO’EO) 

XL=XK0L/YK 

SIOl=Xi<0SIG/X< 

H=2.»PI*P 
THFTA=  0.0 
KK  = 8 

WRITE(6,201)F 

201  FORMATHH  , 2X  . 2H  F=  , F 15  .2  ) 

HPITit6,2C4»XKUL 

204  FOHMATllll  ,2Y,3H<L=,F7.51 
HrITE  (b,20r>»  XKOSIG 

205  FORMAT  IIH  , 2X , GhKSIG 1* ,F 7 .51 
WkITE(b,20f.)- R 

206  FORMATdH  ,2X,3HEi<  = ,F7.5) 

LL  = 4 

NN  = 4 

no  12  1=1, LL 
00  13  J=l,N'i 
WRITE (6,2C2IZXP 

202  FORMAT (IH  ,?X,4H7XP=,F7.5» 

WPirE16,2''3)7YP 

203  FORMATllH  , 2X , 4HZ YP= , F 7 .5 ) 

CO  14  N=1,K< 
rH2=THRrA*pr/i«o . 

TAMT-T  AN  ITur  ) 

DCQ=l.  /SORKl . ♦ZXP**2*2YP**2) 

QJ=  (0  .C,  l.C) 

QEIX=-COSlIH21 
EIT  = 0 . 0 
nEIZ=-51N(TH2) 

C05rp= (ZXP«3IN  irH2)»C0SlTH2M ’DOO 

SNrSIN (TH7I 

CN=COi<TM2l 

liINTP=3aRT  tl.-COSTo»»21 
IF(ZYP)23u,26r ,250 

25u  01  = 1.7iiQRT(7YP*»2  + (SN-ZXP»CN»  ’’Z) 

QAt=ni*lZYf • lQFIX»CN»uEIZ»ON)YFIY» (SN-ZXP»CN» ) 

OA21=  tSN-ZXP»CNl •IOErZ»SN»QEIX*ON) 

QA2=XK» (ZYF»ET  Y-JA21 ) *01/ {H»U01 
ALXl  = 0UO*Ol*  ( SN-ZXP*CH«-SN*ZYP**2> 

ALX2  = -COO*Ot*ZYP*  (C.N*ZXP*SN1  ■ 


30 


ALKT- J Ja'Ol' (iXP* (SN-ZXP*CN»-CN»ZVP»*2) 

ALYl  = ZV''»rM»01 
ALrZ’Dl* ( SN-ZXP*CN» 

ALY3=T 1*ZVP»5N 
ALZ1=- ZXP*COO 
ALZZ=-ZYP»CD1 
ALZJ=  UCO 
GO  n zzn 

?6u  IF(ZYP)261,2o2,261 

262  UAls'^IY 

QA2r-XK»(QEIZ*SN+QEIX*CU)/(H*UO> 

ALX1  = 1 .0 
ALXZrO  .0 
ALX3=0  .0 
ALY1  = 0 .0 
ALY2  = 1 .0 
ALY3-0  .0 
AL71=0.!! 

ALZ2  = 0 .0 
AL73=1 .0 
GO  ro  270 

261  IF(ZXH-TANT>263,264,263 

263  SNAjr ( SN-ZXP»CM / A0S (SN-ZXP»CN) 

UAlr^-I  Y*SNA3 

QA2  = -XK»  «QEIZ»SN+QEIX»CNl*SNAB/lM*UO> 

ALXl=!SNA3/Sa'’T  (l.tZXP**2t 
ALXZ^O  .0 

ALX3=ZXP»SNA3/saRT(l.+ZXP»*Z» 

ALYl^O  .0 
ALY2ISNAB 
ALY3  = 3 .0 

ALZ1=-  ZXP /SiiRT  (1 
ALZ2=0  .A 

ALZ3=  1 ./SaPT(  1 .«-ZXP**2) 

GO  TO  270 

264  QAls'^IY 

QA2  = -XK*  ( 0''IZ»SN  + QErx*CN)/  tW’UO) 

ALX1=1 ,/SORT( 1 .♦TANT**2) 

ALX2=0  .0 

ALX3  = TANr/SQRr  (l.+TANT**2T 
ALY1  = 0 .0 
ALY2=1  .0 
ALY3  = 0 .0 

AUZ1  = -TANT/SQPY11 . ♦TAnT**2» 

ALZ2=J.0 

ALZ3=1./SQPT<1,»TANT»»2) 

270  XKP=H»SQPT(U0»£0*EO> 

P1=CJSTP»XKP»»2-X<»SQRT (XKP»»2-(X<»SINTP )*»2) 
«2=COSTP»XKP»*2*XK*SGRT (XKP»*2-(XK»SINrp) **2) 

RPP=R1/R2 
TPP=1. +RPP 

QCO  = UJ*(XK*»?  -XKP**2)»  <1.-RPPI  'r,0STP/(2.*PI*XK) 
QCl»aJ*(2.»X<*SINTP*SINTP)  *TPP*(XK**2-XKP**2I  /(2.*PI»XKP*'»2) 
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LrAi3]-fr(uu/CD) 

XKc’-X<»rQSTP 

XK2>'  = SU'<T(XK’»»2-«XK»SINTP)»»2I 

CIC2  = -UJ»  rp.-»XK’SINrP*SI.'XTP»|XK»»2-XKP»»2>/<2.*PI*XKP*»2l 

aflll:-XK2»((XK*XK2PI **2) 'Qro 

U*22  = - XK2»XK2  ^*XK»XKP*XK"»OCl 

QA7i--XVrZ»XK2^*X<»XKP»X<P»UC2 

QAi.i.  = X kZP»  (x<Z'X<p»»»2*A<Z»CXK»XKZn»  ’’Z 

QA-  ( jAll*aA22*QA33)/QA<.i* 

QAA11  = QA2»CTA*QA»C J3TP 

»P1  = XK •COSTP-SCPT 1 X<P»»2-C  XK»STNTP»**2I 
RP2  = XK»COSTPf SQKr ( XKP • *2- < XK»S INTP  l•*2l 
RPiRP 1 /PP? 

TP=1. fOP 

aAA22  = -Qj*QAl»  TO» ( XKP* *2-x K* *2 » / (2 .’PI*  <XKZ*XKZP»I 

QAAAX  = CA?»I.A*ETA<  SINTP 

UAXH=aAAll*ALXl+uAA22*ALyi*QAA35»ALZl 
GAYH -n.AAU*ALX2*QAA22»ALY2*QAA33»ALZ2 
QAZH^U^All* ALXo»QAA22»ALY3 ♦aAA33»ALZ3 
QFPX:- CM 

ERYrn. 0 

QIPZ--SN 

QALP=CAXH»0EPX ♦3AYH*eRY»QAZH»GERZ 
UALPn=rOMJGI jALP» 

W3  = '1.3*XL»yL'tXP(  - (XK*XL*SINTP)  »*2) 

QSIiJA=3. * (PI •X<«SIG1»C0STP) *»2»UALP*QALPC*HS 
STr,0A=  PEAL  IQSI  G3A) 

SIGt'f1  = 10  . ‘ALGilO  (SIGOA) 
writ:  t 6,2  0 7)  T HETA,  oIGOd 

207  FOR.IAT (/2X,6HTHETA=,F5.2,10X,17HSIGMA  ZERO  IN  OB=iF2.2) 
THrir  A=  TPETA+l  C . 

14  CONTINUE 
THETA=  0. 0 
ZXP=ZXP*-0.4 
13  CONTINUE 
THFTA=C. 0 
ZXP  = 0 . 0 
ZYP=ZYPf 0.4 
12  CONTINUE 
STOP 
ENO 
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i>  tho  square  rout  of  minus  one  ( vTT, 

IS  frequency, 
is  27rf. 

is  tune  in  seconds. 

is  a unit  vector  in  tlie  direction  of  propagation, 
is  the  angle  of  incidence. 
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Z (X.  y) 
s(x,  y) 


are  unit  vectors  in  tlie  x.  y.  and  z directions,  respectively. 

is  tlie  local  angle  of  incidence. 

is  the  expression  for  the  tilted  plane. 

is  the  expression  for  the  random  slightly  rough  surface. 


U2 


Ei 

k 

O 

« 

r 


1 

a , a. 

-.X 

-y  ’ - z 

e. 

, e.  , e. 

IX 

Hi 

’ ly  ’ 1 

Mo 

k , 

,k 

X ’ 

y 

V.  (kx’k  I’ 


k' 


is  a unit  vector  which  is  nonnal  to  the  surface  Z (x.  y). 

is  a unit  vector  in  the  direction  of  the  receiver. 

is  the  slope  of  the  tilted  plane  in  the  x direction. 

is  the  slope  of  the  tilted  plane  in  the  y direction. 

is  the  incident  electric  field. 

is  Mie  free  space  propagation  constant. 

is  a position  vector. 

are  unit  vectors  in  the  x’.  y',  and  i!  directions,  respectively. 

is  the  unit  polarization  vector  for  the  incident  wave. 

are  the  components  of  gj  in  the  x,  y,  and  z,  directions  respectively. 

is  the  incident  magnetic  field  vector. 

is  the  permeability  of  free  space. 

are  the  Fourier  transform  variables. 

A , (k  , k ),  A , (k  , k ) are  the  amplitudes  of  the  scattered  waves  in 
yl  X y zi  X y 

Fourier  transform  variables  for  a local  horizontally  polarized  incident 
wave. 

is  the  Fresnel  reflection  coefficient  for  a horizontally  polarized  incident 
wave. 

is  the  propagation  constant  in  the  medium  below  the  surface, 
is  the  relative  dielectric  constant  of  the  medium  below  the  surface, 
is  the  permittivity  of  free  space. 

is  the  transmission  coefficient  for  a horizontally  polarized  incident 
wave. 
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A , (k  . k ),  A , (k  , k ),  A , (k  , k ) are  tlie  amplitudes  of  the  scattered  waves  in 
\ ^ X y X y 

Fourier  translorm  variables  for  a local  vertically  polarized  incident 
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is  the  intrinsic  impedance  of  free  space, 
is  the  speed  of  littht  in  free  space. 

is  the  Fresnel  rellection  coefficient  for  a vertically  polarized  incident 
wave. 

is  the  transmission  coefficier.t  for  a vertically  polarized  incident  wave, 
is  the  backscattered  far  field. 

is  the  distance  from  the  origin  of  the  local  coordinate  system  (x'.y'.z'l 

to  the  receiver  in  the  far  field. 

is  the  unit  polarization  vector  for  the  received  wave. 

is  the  urea  on  the  surface  which  is  illuminated. 

is  the  correlation  function  of  the  random  rough  surface. 

is  the  variance  of  the  undulations  of  the  random  rough  surface. 

is  the  normalized  correlation  function  of  the  random  rough  surface. 

is  the  surface  roughness  spectrum  equal  to  the  Fourier  transform  of 

C(T^,Ty). 

is  the  correlation  distance  associated  with  C ). 

is  the  radar  backscatter  coefficient. 
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